ON THE COMPLEXITY OF PUTLNAR'S POSITIVSTELLENSATZ 
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Abstract. Let S = {x g R n \ gi(x) > 0, . . . ,g m (x) > 0} be a basic closed 
semialgebraic set defined by real polynomials g;. Putinar's Positivstellensatz 
says that, under a certain condition stronger than compactness of S, every real 
polynomial / positive on S posesses a representation / = 5Z™=0 a i9i wnere 
go := 1 and each cr; is a sum of squares of polynomials. Such a representation 
is a certificate for the nonnegativity of / on S. We give a bound on the degrees 
of the terms <J{gi in this representation which depends on the description of 
S, the degree of / and a measure of how close / is to having a zero on S. As 
a consequence, we get information about the convergence rate of Lasserre's 
procedure for optimization of a polynomial subject to polynomial constraints. 



1. Introduction 

Always write N := {0, 1,2,...} and R for the sets of nonnegative integers and 
real numbers, respectively. Denote by M[X] the ring of polynomials in n > 1 
indeterminates X := (Xi, . . . , X n ). We use suggestive notation like R[X] := {p 2 \ 

— _ 2 

p £ R[X]} for the set of squares and 2JR[X] for the set of sums of squares of 
polynomials in R[X]. A subset M C R[X] is called a quadratic module if it contains 
1 and it is closed under addition and under multiplication with squares, i.e., 

1 £ M, M + M C M and R[X} 2 M C M. 

A subset T C R[A > ] is called a preordering if it contains all squares in R[A] and it 
is closed under addition and multiplication, i.e., 

R[A] 2 C T, T + T C T and TT C T. 

In other words, the preorderings are exactly the multiplicatively closed quadratic 
modules. 

Throughout the article, we fix m £ N and a tuple g :— (g±, . . . , g m ) of polynomials 
gi £ R[X]. It will be convenient to set go := 1 £ M[X]. The quadratic module M(g) 
generated by g (i.e., the smallest quadratic module containing each gi) is 

m 
i=0 



(1) M(g) = ^^R[A] 2 ffi := \ I ^ e ^R[X]^ 

Using the notation 




-5 <5i 5 m 
9 ■= 9i ■••9m, 



Date: December 15, 2008. 

2000 Mathematics Subject Classification. Primary 11E25, 13J30; Secondary 14P10, 44A60, 
68 W40, 90C22. 

Key words and phrases. Positivstellensatz, complexity, positive polynomial, sum of squares, 
quadratic module, moment problem, optimization of polynomials. 

The first author is supported by National Science Foundation DMS-0456960. 
The second author is supported by the DFG grant "Barrieren" . 

1 



2 



JIAWANG NIE AND MARKUS SCHWEIGHOFER 



the preordering T(g) generated by g can be written as 

( 2 ) T (g)= E £ w 5 k^E R H 2 • 

5e{o,i} m [<5e{o,i} m J 

i.e., T(g) is the quadratic module generated by the 2 m products of gi. It is obvious 
that all polynomials lying in T(g) D M{g) are nonnegative on the set 

S(g) := {x e K" | ^(z) > 0, . . . , ffm (x) > 0}. 

Sets of this form are important in semialgebraic geometry (see [BCRj ) and are called 
basic closed semialgebraic sets. In 1991, Schmudgen [Smnj proved the following 
"Positivstellensatz" (a commonly used German term explained by the analogy with 
Hilbert's Nullstellensatz). 

Theorem 1 (Schmudgen). Suppose the basic closed semialgebraic set S(g) is com- 
pact. Then for every polynomial f G M[X], 

/ > on S(g) / S T{g). 

Under a certain extra property on M{g) which we will define now, this theorem 
remains true with T{g) replaced by its subset M(g). We introduce the notation 

71 

||x|| 2 :=E^ 2 eRW- 

z=l 

Definition 2. A quadratic module M C M.[X] is called archimedean if 

N - \\X\\ 2 e M for some JVeR 

Note that this definition applies also to preorderings since every preordering is a 
quadratic module. As a corollary from Schmudgen's Theorem, we get the following 
well-known characterization of archimedean quadratic modules. 

Corollary 3. For a quadratic module M C K[X], the following are equivalent. 

(i) M is archimedean. 

(ii) There is a polynomial p G M such that S(p) = {p > 0} C M n is compact. 
(Hi) There is a tuple g of polynomials such that S(g) is compact and M contains 

the preordering T(g). 
(iv) For all p G R[X], there is N G N such that N -pE M. 

Proof. Observe that Q ==>• ((uj (pTI| (pv)) =^> (P . All of these implications 
are trivial except (|m)) =>■ ([rv| which follows from Theorem [1] □ 

In particular, we see that S(g) is compact if and only if T(g) is archimedean. 
Unfortunately, S(g) might be compact without M(g) being archimedean (see |PD[ 
Example 6.3.1]). What has to be added to compactness of S(g) in order to ensure 
that M(g) is archimedean has been extensively investigated by Jacobi and Prestel 
|JP[|PD| . Now we can state the Positivstellensatz proved by Putinar [Put] in 1993. 

Theorem 4 (Putinar). Suppose the quadratic module M(g) is archimedean. Then 
for every f G R[X], 

f > on S(g) ==> / G M(g). 
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Both the proofs of Schmiidgen and Putinar use functional analysis and real 
algebraic geometry. They do not give information how to construct a representation 
of / showing that / lies in the preordering (an expression like in ^ involving 2 m 
sums of squares) or the quadratic module (a representation like in |T]) with m + 1 
sums of squares). 

Based on an old theorem of Polya [Pol] , new proofs of both Schmudgen's and 
Putinar's Positivstellensatz have been given in |Swl|ISw3] which are to some extent 
constructive. By carefully analyzing a tame version of |Sw3] and using an effec- 
tive version of Polya's theorem | PRJ , upper bounds on the degrees of the sums of 
squares appearing in Schmudgen's preordering representation have been obtained 
in |Sw2j . The aim of this article is to prove bounds on Putinar's quadratic module 
representation. They will depend on the same data but will be worse than the ones 
known for Schmudgen's theorem. 

Since it will appear in our bound, we will need a convenient measure of the size 
of the coefficients of a polynomial. For a G N n , we introduce the notation 

\a\ := ai + • • • + a n and X a := Xf 1 • • • X%« 

as well as the multinomial coefficient 

\a\\ \a\\ 



a. J a.\\ . . . a n \ 

For a polynomial / = ^ Q a a X a G M[X] with coefficients a a G R, we set 

„,,, \a n \ 
If := max- 



OS 1 ) 

This defines a norm on the real vector space M[X] with convenient properties il- 
lustrated by Proposition If 41 below. For any k G R>o, we now define convex cones 
T(g,k) and M(g,k) in the finite-dimensional vector space R[X] <fe of polynomials 
of degree at most k (i.e., at most [&J) by setting 

T(g,k)=l Yl <^89 S \<J 5 eJ2nX] 2 ,deg(a 5 g s )<k\ C T{g) n R[X]< fc , 
[<5e{o,i}™ J 

M(g,k) = I ^ e ^M[X] 2 ,deg(a 5 3 5 ) < fcj C M{g) nl[I]< fe 

We now recall the previously proved bound for Schmudgen's theorem. 

Theorem 5 ( Sw2 ). For all g defining a basic closed semialgebraic set S(g) which 
is non-empty and contained in the open hypercube (—1,1)", there is some c > 1 
(depending on g) such that for all f £ M.[X] of degree d with 

f* := min{/(x) | x G S(g)} > 0, 

we have 

f€T(g,cd 2 (l+(dwM 



In this article, we will prove the following bound for Putinar's theorem. 
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Theorem 6. For all g defining an archimedean quadratic module M(g) and a set 
7^ S(g) C (—1,1)™, there is some c G R>o (depending on g) such that for all 
f G WL[X] of degree d with 

f* := mm{f{x) | x G S(g)} > 0, 

we have 

feM ^, ceX p^dvwy)). 

In both theorems above, there have been made additional assumptions compared 
to Schmudgen's and Putinar's original results. But these are not very serious 
and have only been made to simplify the statements: For example, if S(g) = 0, 
then — 1 G T(g,k) for some k G N by Schmudgen's theorem. Therefore 4/ = 
(/+l) 2 + (/-l) 2 (-l) e T(g,2d+k) for each / G R[X] of degree d > 0. The other 
hypothesis that S(g) be contained in the open hypercube (—1, l) n is only a matter 
of rescaling by a linear (or afhne linear) transformation on M n . For example, if r > 
is such that S(g) C (— r, r) n , then Theorem [5] remains true with ||/|| replaced by 
||/(rX)||. Here it is important to note that the property that M(g) be archimedean 
is preserved under affine linear coordinate changes. This is clear from Corollary [3] 
Confer also the proof of Proposition [9] below. 

In both Theorem [5] and [6l the bound depends on three parameters: 

• The description g of the basic closed semialgebraic set, 

• the degree d of / and 

• a measure of how close / comes to have a zero on S(g), namely ||/||//*. 
The main difference between the two bounds is the exponential function appearing 
in the degree bound for the quadratic module representation. It is an open research 
problem whether this exponential function can be avoided. It could even be pos- 
sible that the same bound than for Schmudgen's theorem holds also for Putinar's 
theorem. In view of the impact on the convergence rate of Lasserre's optimization 
procedure (see Section [2] below) , this question seems very interesting for applica- 
tions. Whereas the bound for the preordering representation cannot be improved 
significantly (see [Stej ). this seems possible for the quadratic module representation. 

The dependance on the third parameter ||/||//* is consistent with the fact that 
the condition /* > cannot be weakened to /* > in neither Schmudgen's nor 
Putinar's theorem. Under certain conditions (e.g., on the derivatives of /), both 
theorems can however be extended to nonnegative polynomials (see |Sch[ IMr2| ). 
With the partially constructive approach from [Sw4j to representation of nonneg- 
ative polynomials with zeros, one might perhaps in the future gain bounds even 
for the case of nonnegative polynomials which depend however on further data (for 
example the norm of the Hessian at the zeros). 

In special cases, Prestel had already proved the mere existence of a degree bound 
for Putinar's Theorem depending on the three parameters described above (see 
[PD| Section 8.4] and [Pre] ). He used model theory and valuation theory to get 
the existence of such a bound. But the only information about the bound he gets 
(using Godel's theorem on the completeness of first order logic) is that the bound 
is computable. 

In contrast to this, our more constructive approach yields information in what 
way the above bound depends on the two parameters d and ||/||//*- The constant 
c depends on the description g of the semialgebraic set, but no explicit formula is 
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given. For a concretely given g, one could possibly determine a constant c like in 
Theorems [5] and E] by a very (probably too) tedious analysis of the proofs (cf. |Sw2| 
Remark 10]). 

We conclude this introduction by considering the one variable case, i.e., n = 1. 
Scheiderer showed in |Sch[ Corollary 3.4] that, in this case, compactness of S(g) 
implies that M(g) = T(g) (and therefore M(g) is archimedean) . Now the equality 
M(g) = T(g) implies in particular that g 5 G M{g) for all 5 G {0, As an easy 
consequence, we get that Theorem [5] remains valid with T replaced by M in the 
case of univariate polynomials. The bound in Theorem [5] is thus far from being 
sharp in the one variable case. As said above, in the multivariate case it is not 
known if the bound can be improved considerably. 

The rest of the paper is organized as follows. In the next section, we use our 
result to investigate the accuracy of Lasserre's "sums of squares relaxations" for 
optimization of polynomials. In Section [31 we give the proof of Theorem [5] 

2. Convergence rate of Lasserre's procedure 

Consider the problem to compute (by a numerical procedure, i.e., up to some 
prescribable error) the minimum 

(3) f* ■■= mm{f(x) | x G S(g)} 

of a polynomial / G M[X] on a non-empty basic closed semialgebraic set S(g). 
In other words, you want to minimize a polynomial under polynomial inequality 
constraints. When all the polynomials involved are linear, i.e., of degree < 1, this 
is a linear optimization problem (a linear program) and there are very efficient 
algorithms to solve this problem. For general polynomials this problem gets very 
hard. It is therefore a common approach to solve a much easier related problem, a 
so called relaxation, namely to compute for k G N, 

(4) := sup{a G E | / - a G M(g, k)} etU {-oo} 

which is clearly a lower bound of /*. The problem of finding f% can be written as a 
semidefinite program whose size gets bigger when k grows (see the references below) . 
Semidefinite programming is a well-known generalization of linear programming for 
which very efficient algorithms exist (see for example |Todj ) . One can now solve 
a sequence of larger and larger semidefinite programs in order to get tighter and 
tighter lower bounds for /*. Lasserre |Lasj was the first to interpret Putinar's 
theorem as a convergence result. 

Indeed, it is easy to see that Putinar's theorem just says that the ascending 
sequence (f^)kefi converges to /* under the condition that M(g) be archimedean. 
In this section, we will interpret our bound for Putinar's Positivstellensatz as a 
result about the speed of convergence of this sequence. 

For an introduction to the interplay of semidefinite programming, sums of squares, 
optimization of polynomials and results about positive polynomials, we refer to 
[Lasi IMrl, Swl] (with special regard to Putinar's Positivstellensatz) and fJL, DNP, 
INDSl IPS] . There are several software tools which translate the problem of com- 
puting into a semidefinite program and call a semidefinite programming solver. 
See pllKKWllLofllSoS] . 

The following technical lemma will also be needed in Section [3l 
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Lemma 7. For any polynomial f £ M[X] of degree d > 1 and all x £ [— 1, 1]", 

|/(ar)|<2dn d ||/||. 

Proof. Writing / = J2 a a a ('"')^" Q ( a « e K), we have ||/|| = max a |a Q | and 



< 



5>«i( l a>>' < " 



for all x S [—1, 1]™. Using that \a a \ < \\f\\ and |xj| < 1, the multinomial identity 
now shows that \f(x)\ < ||/||£fc =0 n fc < (d+l)n d \\f\\ < 2dn d ||/||. □ 

Now we are ready to prove the main theorem of this section. 

Theorem 8. For all polynomials g defining an archimedean quadratic module M(g) 
and a set ^ S(g) C (—1,1)", there is some c > (depending on g) such that 
for all f £ WL[X] of degree d with minimum f* on S and for all integers k > 
cexp((2d 2 n d ) c ), we have 

log 7 



and hence 

osr-Jis"*"" 1 



log; 



where f£ is defined as in Q). 

Proof. Given g, we choose c > like in Theorem [5] Now let / £ M[X] be of degree 
d with minimum /* on S and 

(5) k > cexp((2dV) G ) 

be an integer. The case d = is trivial. We assume therefore d > 1. Note that 
> c and hence log(fc/c) > 0. Setting 

_ 6d 3 n 2d \\f\\ 



log 7 



(G) 



all we have to prove is h :=/—/* + a € M(#, fc) because the second claim follows 
from this. By our choice of c and the observation deg h = deg / = d, it is enough 
to show that 

cexpffdVMY) <fc, 

or equivalently 



dV||ft|| < aylog^ = 6d 3 n 2d ||/||. 
Observing that \\h\\ < \\f\\ + |/*| + a, it suffices to show that 

11/11 + I/I + a<6dn d ||/||. 
Lemma [7] tells us that |/*| < 2dn d ||/|| and we are thus reduced to verify that 

a<(4dn d ~ 1)||/|| 
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which is by (J6j) equivalent to 



6d 3 n 2d < (Adn d - l)tflog*. 
By ©, it is finally enough to check that 6d 3 n 2d < (Adn d ~ l)(2d 2 n d ). □ 

As already said in the introduction, the hypothesis that S(g) is contained in the 
open unit hypercube is just a technicality to avoid that the bound gets even more 
complicated. In fact, if one does not insist on all the information given in Theorem 
[51 one gets a corollary which is easy to remember and still gives the most important 
part of information. 

Corollary 9. Suppose M(g) is archimedean, S(g) ^ and f G WL[X]. There is 

• a constant c > depending only on g and 

• a constant d > depending on g and f 

such that for f* and f£ as defined in ^j) and O), 

d 

< /* - fk < -- for all large fc G N. 

V lo g| 

Proof. Without loss of generality, assume / ^ 0. Set d := deg/. Since M(g) 
is archimedean, S(g) is compact. We can hence choose a rescaling factor r > 
depending only on g such that S(g(rX)) C (—1, 1)™. Here g{rX) denotes the tuple 
of rescaled polynomials gi(rX). Now Theorem[8japplied to g(rX) instead of g yields 
c > that will together with d := 6d 3 n 2d ||/(rX)|| have the desired properties by 
simple scaling arguments. □ 

Remark 10. The bound on the difference /* — f£ presented in this section is much 
worse than the corresponding one presented in |Sw2[ Section 2] which is based on 
preordering representations (i.e., where would be defined using T(g) instead of 
M{gj). This raises the question whether it is after all not such a bad thing to use 
preordering (instead of quadratic module) representations for optimization though 
they involve the 2 m products g s letting the semidefinite programs get huge when 
m is not small. However, it is not known if Theorem [8J holds perhaps even with the 
bound from [Sw2[ Theorem 4] . Compare also [Sw2[ Remark 5] . 



3. The proof 

In this section, we give the proof of Theorem [6l The three main ingredients are 

• the bound for Schmiidgen's theorem presented in Theorem [5] above, 

• ideas from the (to some extent constructive) proof of Putinar's theorem in 

Sw3, Section 2] and 

• the Lojasiewicz inequality from semialgebraic geometry. 

We start with some simple facts from calculus. 
Lemma 11. IfO ^ / S M[X] has degree d, then 

\f(x)-f(y)\ < Hz-iHIdV-^ll/ll 

for all x, y G [-1, 1]™. 
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Proof. Denoting by Df the derivative of /, by the mean value theorem, it is enough 
to show that 



(7) 



\Df(x)(e)\<<Pn d - 1 tfi\\f\\ 



for all x E [—1,1]" and e G M" with ||e|| = 1. A small computation (compare the 
proof of Lemma [7]) shows that 



df(x) 



dxi 



d 

E 

fe=i 



< n/ii YM\ Xl \ + ••• + ! 



from which we conclude for all x G 
^ df{x) 



\Df(x)(e)\ = 



i=i 



dxi 



< 



fc_1 < WfWY^ 1 < \\f\\d 2 n d -\ 

k=l 

1, 1]" and e G E™ with ||ej| = 1, 
df(x) 



E 



dxi 



Because for a vector e on the unit sphere in 
this implies ((7]). 

Remark 12. For all k G N and y G [0, 1], 

(2/ 



•le^H/llrfV-^M. 

2=1 

Ei=i l e *l can reach at most 



'n, 
□ 



i) 2fe y< 



l 



2k + 1 



The next lemma is a version of |Sw3i Lemma 2.3] caring about complexity issues. 
In |Sw3[ Lemma 2.3], it is shown that, if C C W 1 is any compact set, gi < 1 on C 
for all i and / S K[X] is a polynomial with / > on S(g), then there exists A > 
such that for all sufficiently large k S N, 



(8) 



Z-A^te-l) 2 ^ >0 on C. 



i=l 



The idea is that, if you want to show that / £ M(g), you first subtract another 
polynomial from / which lies obviously in M(g) such that the difference can be 
proved to lie in M(g) as well. This other polynomial must necessarily be nonncg- 
ative on S(g) but it should take on only very small values on S(g) so that the 
difference is still positive on S(g). On the region where you are outside and not too 
far away from S(g), the polynomial you subtract should take large negative values 
so that the difference gets positive on this region outside of S(g) (where / itself 
might be negative). The hope is that the difference satisfies an improved positivity 
condition which will help us to show that it lies in M{g). To understand the lemma, 
it is helpful to observe that the pointwise limit for k — > oo of this difference, which 
is the left hand side of (jTTJ) , is / on S(g) and oo outside of S(g). 

Lemma 13. For all g such that S := S(g) n [-1, 1]" ^ and gi < 1 on [-1, 1]", 
there are Cq,C\,C2 > with the following property: 

For all polynomials f G R[X] of degree d with minimum f* > on S , if we set 



(9) L := d 

and if k G N satisfies 
(10) 



2„d-l Wf\ 



\:= Cl d z n a -'\\f\\L c * 



2fc + l >c (l + L Co ), 
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then the inequality 

(11) /-A]>>-l) 2fe ffi >^- 
holds on [— 1, 1]". 

Proof. By the Lojasiewicz inequality for semialgebraic functions (Corollary 2.6.7 in 
BCR]), we can choose 02,03 > such that 

(12) dist(x, S) C2 < -ca mm{gi(x), ...,g m (x),0} 

for all x € [—1, 1]™ where dist(x, S) denotes the distance of x to S. Set 

(13) c 4 :=c 3 (4n) C2 , 

(14) ci := 4nc 4 

and choose Co € N big enough to guarantee that 

(15) c (l + r c °) > 2(m- l)c 4 r C2 and 

(16) c (l + r Co ) > 4mcir C2+1 

for all r > 0. Now suppose / £ K.[X] is of degree d with minimum /* > on S and 
consider the set 

a-.= |xe[-i,in/(x)<^r 

By Lemma [IT] we get for all x £ A and y E S 

f — < f(y) - f{x) <\\x- ylldV-^H/H < Has — y\\d 2 n d \\f\\ 
Since this is valid for arbitrary y £ 5, it holds that 

<dist(^S) 

for all x £ A. We combine this now with (| 1 2[) and get 



min{.gi(a;), . . . ,g m {x)} < 



c 3 V4d 2 ?i d ||/|| 



for x € A. We have omitted the argument in the minimum which is here redundant 
because of A n S = 0. By setting 

(17) S := — L_ > 0, 

where we define L like in |9]), and having a look at (113|) . we can rewrite this as 

(18) min{5i(a;), . . . ,g m {x)} < -5. 
Define A and k like in ((9J and (fT0|) . For later use, we note 

(19) \ = Cl L C2+1 f*. 
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X5 

■ f + Y 


> 


r 


s 


> 


m — 1 


2 


2k + 1 


/* 


> 


Am 


4 


2fc + 1 



We claim now that 
(20) 

(21) " :- 7 <"" 1 

(22) 

Let us prove these claims. If we choose in Lemma 111 I for y a minimizer of / on 
S, we obtain 

\f{x) -f*\< diam([-l, l] n )dPn d ~ l Vn||/|| - 2^d 2 n d - 1 ^\\f\\ = 2d 2 n d \\f\\ 

for all x e [—1, 1]™, noting that the diameter of [—1. 1]™ is 2^/n. In particular, we 
observe 

/ > ./* - 2rf 2 ?1 d ||/|| > f — - 2dV||/|| on [-1, 1]". 



Together with the equation 



Y = 2dV||/||, 



which is clear from ©, (O and (H}, this yields ([20]). Using (fT0| , JT5|) and p7| . 
we see that 

(2fc + 1)5 > co(l + L C0 )<5 > 2(m - l)c 4 L C2 (5 = 2(m - 1) 
which is nothing else than (|2T|) . Finally, we exploit pH|) . and (fT^)) , to see that 
(2* + 1)/* > cq(1 + L C0 )f* > 4m Cl L C2+1 f* = 4mA, 

i.e., holds. 

Now ([20]). (|2T|) and (J22[) will enable us to show our claim ([IT]). Ifi6i, then in 
the sum 

m 

(23) E^*^) - l ? k 9i{*) 

i=l 

at most m — 1 summands are nonnegative. By Remark 1121 these nonnegative 
summands add up to at most (m — 1)/ (2k + 1). At least one summand is negative, 
even < —S by (fTS)) . All in all, if we evaluate the left hand side of our claim (fTTj) in 
a point x d A, then it is 

„, , in — 1 , . „,, XS . f 8 m — 1 \ f* 

> fix) - X h XS > fix) H hA > — . 

- Jy ' 2k + 1 ~ J y ; 2 \2 2k + 11 ~ 2 



>^T b y HSJ >0 by ( f2Tt 

When we evaluate it in a point x e [— 1, 1]™ \ A, all summands of the sum ([23]) 
might happen to be nonnegative. Again by Remark 1121 they add up to at most 
m/(2k + 1). But at the same time, the definition of A gives us a good lower bound 
on f(x) so that the result is 

3 „ , m /* /* Am /* 
> _/* _ \ > __ + L > __. 

~ 4^ 2fc + l~ 2 4 2k + 1 ~ 2 

>0 by 1(22} 

□ 
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Proposition 14. Ifp,q £ M[X] are both homogeneous (i.e., all of their respective 
monomials have the same degree), then \\pq\\ < \\p\\ \\q\\ . For arbitrary s £ N and 
polynomials ^ pi, . . . ,p s £ 1&[X], we have 

lbi-"P.|| < (l + degpi)---(l + degp s )||pi||---||p s ||. 

Proof. The statement for homogeneous p and q can be found in |Sw2| Lemma 8] . 
The second claim follows from this by writing each pi as a sum pi — ^2 k Pik of 
homogeneous degree k polynomials pik ■ Multiply the pi by distributing out all such 
sums and apply the triangle inequality to the sum which arises in this way. Then 
use 

llpifci ■ • -psk s \\ < Ibifeill ■ ■ • WpskJ < lbi|| ■ • • lb a ||. 

Now factor out ||pi|| • ■ • \\p s \\ an d recombine the terms of the sum which now are all 
constant 1. □ 

Lemma 15. For all c\, Ca, C3 > 0, there is c > such that 

c\ exp(c2r C3 ) < cexp(r c ) for all r > 0. 

Proof. Choose any c > ciexp(c22 C3 ) such that C3 < c/2 and C2 < 2 C / 2 . Then for 
re [0,2], 

ci exp(c2r C3 ) < ciexp(c22 C3 ) < c < cexp(r c ) 
and for r > 2 (observing that c\ < c), 

Ciexp(c 2 r C3 ) < cexp(2 c/2 r c/2 ) < cexp(r c ). 

□ 

We resume the discussion before Lemma [13] With regard to (fTT]l , we can for 
the moment concentrate on polynomials positive on the hypercube [—1, 1]™. If 
this hypercube could be described by a single polynomial inequality, i.e., if we had 
[— 1, 1]™ = S(p) for some p £ K[X], then the idea would be to apply the bound for 
Schmiidgen's Positivstellensatz now. The clue is here that p is a single polynomial 
and hence preordering and quadratic module representations are the same, i.e., 
T(p) = M(p). The following lemma works around the fact that [— 1, l] n = S(p) can 
only happen when n = 1. We round the edges of the hypercube. 

Lemma 16. Let S C (-1, 1)" be compact. Then 1 - \ - (Xf d + ... X^ d ) > on 
S for all sufficiently large d £ N. 



Proof. Consider for each KdeN the set 

A d := L £ S I x\ d + ■ ■ ■ + x 2 n d > 1 



This gives a decreasing sequence A\ D A2 2 ^3 3 • • ■ of compact sets whose 
intersection C\°^ =1 Ad is empty by calculus. By compactness, a finite subintersection 
is empty, i.e., A d = for all large d £ N. □ 

Note that in the proof of Putinar's theorem in [Sw3, Section 2] where we were 
not interested in complexity, a different approach has been taken. Condition ([5]) 
has been established for a polyhedron C which is even bigger than the hypercube, 
so big that preordering representations certifying nonnegativity on C can be turned 
into quadratic module representations certifying nonnegativity on the hypercube. 
The advantage was that we could use Polya's theorem [Pol which is much more el- 
ementary than Schmiidgen's theorem. Despite the existence of the effective version 
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[PRj of that theorem of Polya, it seems that establishing positivity on such a big 
polyhedron C is too expensive from the complexity point of view. Though it is not 
so nice, we therefore work here with a rounded hypercube and Theorem [5] instead. 
We finally attack the proof of Theorem [6j 

Proof of Theorem^ By a simple scaling argument, we may assume that ||<7j|| < 1 
and Qi < 1 on [— 1, 1] for all i. According to Lemma ITBI we can choose d £ N such 
that 

P ■= 1 - j- - (X! d + ■■■+ Xl d ) > on S(g). 
do 

By Putinar's Theorem [H we have p £ M(g) and therefore 

(24) p£M(g,d 1 ) 
for some d\ £ N. Choose 4 such that 

(25) l + deggi<d 2 for alii e {1, m}. 

Now we choose Co , c\ , C2 like in Lemma [T3J define L and A like in @ and choose 
the smallest k £ N satisfying (fT0|) . Then 

(26) 2fc + l <c (l + L C0 ) + 2. 

Let C3 > 1 denote the constant existing by Theorem [5] (which is there called c and 
gives the bound for preordering representations of polynomials positive on S(g)). 
Using Lemma [T51 it is easy to see that we can choose C4, C5, cq, C7, c > satisfying 

(27) c 3 2 C3 r 2+2c3 n C3r < c 4 (exp(c 4 r)) 

(28) 2r + 2 Cl r C2+1 4 (1+rCn)+1 < c 5 exp(r C5 ) 

(29) c 4 exp(2c 4 d2r(l + r c ° +3)) < c 6 exp(r Ce ) 

(30) cfca exp(c 3 r C5 + r c «) < c 7 exp(r C7 ) 

(31) c 7 exp(r C7 ) + d 1 < cexp(r c ) 

for all r > 0. Now let / £ M.[X] be a polynomial of degree d > 1 with 

/* := min{/(x) | x £ S(g)} > 0. 
We are going to apply Theorem [5] to 



/-A^-l) 2 ^ 



i=l 



By Lemma [T31 (|lip holds for this polynomial, in particular 
(32) h* := imn{h(x) \ x £ S(p)} > f -. 

By Proposition [T4l and the definition of c?2 in (|2"5"1) . 



(33) \\h\\ < \\f\\ + Xdf +1 

(34) deg h < max{d, (2k + l)d 2 , 1} =: d h . 

By Theorem [5] (respectively the above choice of C3 > 1), we get 

(35) h £ T(p, k h ) where := c 3 d 2 h (l + d\n dh ^ 
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Note that \\h\\/h* > 1 since < h* < h(0) < \\h\\. We use this to simplify the 
degree bound in (I35|) . Obviously 



(36) k h <c 3 di[2d 2 h n d ^ 



< c 3 2 C3 d 2 h +2c3 n C3dh 
by choice of C4 in (|27|) . Moreover, we have 



< c 4 exp(c 4 c4) 



h* 



(37) 



< 



\f\\+\df +1 ) 



l/ll 



j2fc + l rc 2 + l 



2cidf +1 L c 



ft,* " " 1 ' f* 

<2L + 2 Cl df +1 L C2+1 = 2L + 2c 1 L C2+1 d c 2 o{1+L ' :0)+1 < c 5 exp(L Cs ) 

by H32J), (|2BJ), (HSJ) and by the choice of c 5 in 

4 < d(2k + 2)d 2 

< d(c (l + L ca ) + 3)d 2 



It follows that 

(by m) 

(by ©) 



< 2d 2 cfn 



2„rf 



l/ll 



2dn d \\f\ 



ic {l + L c °) + 3) 



<2d 2 d 2 n d M(c (l + L Co ) + 3) 
/ 

<2d 2 nL(c (l + (ni) Co +3)) 
and therefore 

(38) C4exp(c4d/i) < C6 exp((nZ,) C6 ) 

for the constant cq chosen in (f2"5| . We now get 



(by Lemma [7]) 
(by ©) 



fc/i < C4 exp(c4<f/, 



(by ©) 
(by and 



(by choice of c 7 in (l30l) 'l 



< c 6 exp((nL) C6 )(c 5 exp(L C5 )) C3 
= c^c 6 exp(c 3 (nL) C5 + (nL) c «) 

< c 7 exp((nL) C7 ) 

Combining this with (133)) and (|24p. i.e., 

ft, S T(p, c 7 exp((ni) C7 )) and p £ M(g, di), 
yields (by composing corresponding representations) 

h £ M(g,cexp((nL) c )) 
according to the choice of c in (f3"Tj) . Finally, we have that 



f = h + \J2(9> - l ? k 9i G M(ff, cexp((nL) c )) 



since 



deg(( 5l - l) 2k g t ) < d h <k h < c 7 exp((nL) C7 ) < cexp((nL) c ) 
by choice of d 2 in (|25j) . dft in (f34)) . fe/, in (f35|) and c in (f3Tj) . 



□ 
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